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ABSTRACT 

We present a multi- wavelength study of GRB 081008, at redshift 1.967, by Swift, ROTSE-III 
and GROND. Compared to other Swift GRBs, GRB 081008 has a typical gamma-ray isotropic 
equivalent energy output (~ 10 53 erg) during the prompt phase, and displayed two temporally 
separated clusters of pulses. The early X-ray emission seen by the Swift/XRT was dominated 
by the softening tail of the prompt emission, producing multiple flares during and after the 
Swift/BAT detections. Optical observations that started shortly after the first active phase of 
gamma-ray emission showed two consecutive peaks. We interpret the first optical peak as the 
onset of the afterglow associated with the early burst activities. A second optical peak, coincident 
with the later gamma-ray pulses, imposes a small modification to the otherwise smooth lightcurve 
and thus suggests a minimal contribution from a probable internal component. We suggest the 
early optical variability may be from continuous energy injection into the forward shock front by 
later shells producing the second epoch of burst activities. These early observations thus provide 
a potential probe for the transition from prompt to the afterglow phase. The later lightcurve of 
GRB 081008 displays a smooth steepening in all optical bands and X-ray. The temporal break 
is consistent with being achromatic at the observed wavelengths. Our broad energy coverage 
shortly after the break constrains a spectral break within optical. However, the evolution of the 
break frequency is not observed. We discuss the plausible interpretations of this behavior. 

Subject headings: gamma-ray burst: individual (GRB 081008) 
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Introduction 



The fast and precise localizations of G amma 



ray B ursts (GRBs) provided by Swift/BAT (jGehrels et al 
2004) in hard X-rays have enabled follow-up ob- 
servations in longer wavelengths during the burst 
or shortly after its cessation for an increasing 
number of events. These observations, both 
from ground and from the narrow-field instru- 
ments onboard Swift, have revealed features that 
were not observed in the w e- Swift era. In X- 
rays, a canonical lightcurve (Nouse k et al. 2006; 
IZhang et all2006HO'Brien et al.ll2006l ) starts with 
a rapidly decaying phase that is too steep to be 
explained by the standard external blast wave 
model. This phase is often followed by a slowly 
decaying period, shallower than expected from 
the same external shock model. Only after this, 
the lightcurves show the normal decay and some- 
times a l ate break, often interp r eted as a jet break 
see e.g. IWillmgale et ail 120071 iLiang et al.l l200l 



Racusin et all l2009t lEvans et all 12009^ as ob- 
served for the pre-Swift bursts. Another common 
but unexpected phenomena are the X-ray fl ares 



( Chincarini et al 



pnenon 



Falcone et al] 120071 ) su 



perimposed on the smoothly decaying segments 
observed for about half of the bursts. Their short 
time-scale variability is hard to produce externally 
at a large distance from the progenitor. All these 
features, perhaps explained by the delayed large 
angle internal emission ( e.g. Kumar fc Panaitescu 

" 20071). 



2000; ILiang et all 2006; Z hang et al 



tinuous injection of energy 
19981: ISari 



con- 



Rees fc Meszaros 



Meszarosl 120001 IZhang et alT T2006: 
20071 ). delayed internal shocks or 



Liang et alJ 

central enginTTcti vitv lohincarini et al 



2007 



Falcone et ail 120071 iKruhler et all 120091 and ref- 



erences therein), suggest associations with central 
engine properties. It is likely that the central 
engine is active for much longer than previously 
believed. 

At optical wavelengths, an initial rising optical 
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counte rpart is detected for a significant fraction of 
bursts dMolinari et al. 2007; iPanaitescu fc Vestrandl 
2008 1 ; lOates et al.ll2009j:lKlotz et al.l2009l:lRvkoff et alJ 
20091: IKruhler et alJ l2009t iGreiner et all 120091 
Guidorzi et alJ 2009T ). The temporal behaviors 



generally agree with predictions of the fireball for- 
ward shock model ( Sari fc Piranl 19991 : Meszaros 
2_006j) before the onset of the shock deceleration. 
Alternatively, a rising afterglow can be attributed 
to emission in a collimated afterglow, viewed from 
outside the initial jet opening angle, when the 
shock decel erates and the relativistic beam ing an- 
gle widens ( Panaitescu fc Vestrand 20081) . A va- 
riety of lightcurve shapes, including ones with a 
wide peak or plateau phase, can be accommo- 
dated in this model by combination of the jet 
geometry and the viewing angle. Such signature 
of early external shock emission is not observed 
in X-ray, because the accompanying X-ray emis- 
sion is likely obscured by the dominating internal 
emission. Optical observations thus provide im- 
portant clues about the onset of the afterglow 
because emission from the external shock appears 
to dominate the optical emission from early times. 

During the burst, optical flares apparently cor- 
related with the hard X- ray emission are ob- 
served for several events (IVestrand et al.l 120051 
20061 IPage et alJ l2007t lYuan et alJl2008t l. High 



time-resolution data plays a key role in studying 
the real correlation between the low and high en- 
ergy emission. Such observations in optical, how- 
ever, are only feasible for either a once-in-a-decade 
event, like the "naked eye" burst GRB 080319B 
( Racusin et al.l [2008b ). which was bright enough 
to be detected by a surveying very-wide-field op- 
tical instrument, or a bright burst lasting long 
enough in a usual follow-up scenario (including 
when the Swift/BAT is triggered on a pre-cursor) . 
For a number of other events, an optical flare unre- 
lated to the high e nergy emission is detected (e.g. 
Akerlof et al1ll999l ) or a monotonic decaying tran- 



sient is observed at a flux level consistent with 
back-extrapolation from later observations (e.g. 



Rvkoff et al J 12005). These provide evidence that, 



at least in some cases, emission from the external 
shock dominates the prompt optical observations. 
Such diverse behaviors suggest that both internal 
and external shock emission may contribute to 
the prompt optical detections but their relative 
strengths vary from event to event. 
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GRB 081008, at a redshift of 1.967, is a typi- 
cal long-lasting burst detected by Swift/BAT that 
provides a rare opportunity to study the optical 
characteristics during the prompt phase. ROTSE- 
IIIc started imaging only 42 s after the burst trig- 
ger and before the second epoch of major 7-ray 
emission. An initially rising optical transient was 
observed, followed by two peaks, the latter coinci- 
dent with a gamma-ray peak. Swift slewed imme- 
diately, and the XRT and the UVOT started ob- 
serving during the second optical peak, providing 
well-sampled data at early times and followed the 
event until it dropped below the detection thresh- 
old (see FigureQ}. At about 13 ks after the burst, 
GROND started observing the afterglow in its 7 
optical/IR channels. The excellent energy cover- 
age from IR to X-ray allows us to model the spec- 
tral energy distribution with an intrinsic broken 
power-law, although the evolution of the spectral 
shape is not restricted. The overall behavior of the 
observed afterglow is consistent with being achro- 
matic, while there is some hint of a slightly steeper 
temporal decay in the X-rays after the break. 
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Fig. 1. — Multi- wavelength lightcurve of GRB 
081008. The UVOT flux densities are scaled to the 
u-band using the normalizations determined from 
the simultaneous temporal fit. The ROTSE-III 
unfiltered data are scaled to the GROND r'-band 
using the normalization from the same temporal 
fit. The exponential t o power-law decay models 
(jWillingale et al.l 120071) for the afterglow compo- 
nents are over-plotted in dashed lines. Selected 
epochs or time intervals for spectral analysis are 
labeled with corresponding indices in Figure [5] 



In the following sections, we first summarize the 
Swift, ROTSE-III and GROND observations in § 
2. We then present the multi- wavelength spectral 
and temporal analysis at selected epochs during 
the prompt and afterglow phases in § 3. We next 
discuss the interpretations of the observations in 
§ 4 and finally conclude in § 5. Throughout the 
paper, we adopt the convention that the flux den- 



sity can be described as /„ oc t 



where a 



and P are the corresponding temporal and spec- 
tral power-law indices. Note that a negative a 
represents a rising lightcurve while a positive /3, 
equal to the photon index minus one, corresponds 
to a spectrum that falls with increasing energy. 
We assume a standard cosmology model with the 
Hubble parameter Ho=70 km s _1 Mpc -1 and the 
density parameters f2 m =0.3, Oa=0.7. All quoted 
errors are 1-sigma (68% confidence), unless other- 
wise stated. 

2. Observations and Data Reduction 
2.1. BAT 

Swift BAT was triggered by G RB 081008 on 



19:58: 09.4 UT (Trigger 331093, iRacusin et al 
2008a). We refer to this time as the trigger 
time, Ttriggcr- The BAT data were analyzed us- 
ing the standard analysis software distributed 
within FTOOLS, version 6.5.1. The command 
batgrbproduct was first run to provide the basic 
set of products. This analysis found the emission 
to be detected from T tr ig g or-65.2 s and we there- 
fore adopt this time as the onset of the burst To- 
Such a shift of the reference time only affects our 
analysis of the prompt and the very early after- 
glow, but we use it throughout the paper to be 
consistent. 

Figure[T]shows the BAT flux density (at 71 keV) 
interpolated using the mean spectral index. The 
counts are binned with a minimum signal-to- 
noise ratio of 6 and maximum bin size of 20 s. 
The lightcurve displays a weak component de- 
caying from To, two bright structured peaks at 
Tq+~65 s and To+~175 s and a fainter peak at 
To+~125 s. The emission drops below the BAT 
detection threshold at T o +-270 s. The T90 (15- 
350 keV) of the burst is measured to be 185 ±39 s. 
Spectral hardening and subsequent softening are 
observed during the first broad peak at T o +~70 s. 
The average spectrum of the burst (from Tq to 
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To+263.3 s) is well fit by a power-law with expo- 
nential cutoff model. The best fit photon index 
and cutoff energy are L26(-0.32,+0.24) and 117(- 
50, +83) keV (with 90% confidence). At a redshift 
of 1.9 67 (|Cucchiara et al.ll2008bl; iD'Avanzo et~ail 
2008), the isotropic equivalent energy release 
(0.1 keV-10 MeV) is -6.3xl0 52 erg, given a total 
fluence estimate of ~6.7xl0~ 6 erg cm~ 2 . 

2.2. XRT 

Swift slewed to this burst immediately and the 
XRT began observing in window timing (WT) 
mode at To+152.3 s. A bright uncatalogued X- 
ray counterpart was identified and its position dis- 
tributed via a GCN alert notice. As the count rate 
dropped, XRT switched to photon counting (PC) 
mode at T +476 s. 

The XRT lightcurve was obtaine d from the 
Swift /XRT GRB lightcurve repository (|Evans et al 



20071 ). The spectra were extracted from the XRT 



te am repository, from the new products outlined 



(|2009h . The X-ray flux densities at 
2.77 keV were calculated using mean photon in- 
dices for WT and PC mode accordingly. This cho- 
sen frequency is the weighted mean of the X-ray 
emission assuming /3=1, a typical spectral index 
observed, which is also similar to the afterglow 
spectral index of GRB 081008. The interpolations 
are thus not very sensitive to small changes in the 
photon index and the data from the two modes 
join smoothly. Three peaks are seen in the XRT 
lightcurve in Figure [TJ The first two have complex 
structures and the first peak is at a similar time 
as the last BAT detected pulse. The third peak is 
a flare superimposed on the steep decay after the 
second peak. After that, the lightcurve smoothly 
declined and steepened at To+~15 ks. 

As observed for the first bright BAT peak, 
the XRT hardness ratio (1.5-10 keV/0.3-1.5 keV) 
curve shows a similar structure as the early 
lightcurve 

the flares (see iBurrows et al 

120061 : iRomano et al.l I2006L for other examples of 
such behavior). In general, there is also a trend of 
softening over time. This behavior, together with 
the rapid variation, indicates that the early X- 
ray detections are dominated by continuing burst 
activity, e.g. emission produced in the internal 
shocks. Detailed spectral analysis will be pre- 
sented in § 3. 



dis playing visible hardening during 
1 20051 iFalcone et al.l 



2.3. UVOT 

UVOT began observations of the field at 
To+142.7 s. The observing sequence started with 
a 10 s v band settling exposure and a 150 s white 
(160-650 nm) finding chart exposure. These were 
followed by a 250 s u band exposure, after which 
the filter wheel rotated through all filters, taking 
short 20 s exposures. One more 150 s white band 
finding chart observation was taken at To+792 s. 
After 5 ks, longer exposures were obtained in all 
filters. 

The optical counterpart was detected in white, 
v, 6, u and uvwl bands, but not in uvm2 and 
uvw2 bands. This is consistent with the Lyman-a 
break at the measured redshift. In Figure [TJ all 
UVOT observed fluxes are scaled to u band using 
the normalizations determined from the simulta- 
neous temporal fit described in §3.4. Initial anal- 
ysis of the late time white filter data shows the 
decay o f the afterglow becomes v ery shallow after 
100 ks. ICucchiara et al. ( 2008a ) pointed out the 
presence of an extended source 2 arcsec away from 
the afterglow, which is within the 3 arcsec extrac- 
tion radii used for UVOT. To measure the effect 
of this probable host galaxy, target-of-opportunity 
observations in the white filter were requested on 
Nov. 13th. After removing the host contribution 
from the white band photometry, no detections af- 
ter 90 ks remain above 3 sigma. Such a correction 
is not attempted in other bands, given the negligi- 
ble contribution from the underlying host galaxy 
during these earlier observations. 

2.4. ROTSE-III 



ROTSE-IIIc (|Akerlof et al.l 120031 ). a 45 cm 



robotic optical telescope located at the H.E.S.S. 
site in Namibia, responded to GRB 081008 au- 
tomatically upon receiving the GCN notice. The 
first 5-seconds exposure started 107.1 s after T 
(8.3 s after receiving the alert). An initially bright- 
ening op tical counterp art was detected and re- 
ported in lRvkoffl (|2008h . 

ROTSE-IIIc took a total of 10x5 s, 10x20 s 
and 144x60 s exposures before the target eleva- 
tion dropped below the observing limit. The first 
ten 5-seconds exposures were taken in sub-frame 
mode with 3 seconds readout gap inbetween, in- 
stead of the 6 seconds full-frame readout time for 
later images. We took the first 70 single images 
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and co- added the later images into sums of 5 to 10 
depending on the image quality. The ROTSE-IIIc 
images were bias-subtracted and flat-fielded by 
the automated pipeline. Initial object detections 
were performed by SExtractor ( Bertin fc Arnoutsl 
1996). The images were then processed with our 
custom RPHOT photometry p rogram based on 
the DAOPHOT (IStetsonl 119871) PSF-fi tting pho- 
tometry package (jOuimbv et al. 20061 ) . A refer- 
ence image was constructed from images taken on 
October 18th when the OT had faded below our 
detection limit. 

Due to ROTSE-III's large pixel scale (3.4 arc- 
sec), the optical counterpart is slightly blended 
with four surrounding objects in the ROTSE- 
IIIc images. To remove the possible contamina- 
tion, we tried several different methods. First, we 
used the RPHOT functionality to subtract four 
PSF-scaled point sources around the OT. The 
lightcurve did appear to be steeper than obtained 
without subtraction. The two objects north of 
the OT are blended sources themselves but not 
resolved by ROTSE-IIIc. They may not be com- 
pletely removed using the point source approxi- 
mation, but the residual should be negligible com- 
pared to the bright OT. To confirm this, we at- 
tempted an alternative su btraction method de- 
scribed in lYuan fc Akerlol ( 2008 ) . Each burst im- 
age was cross-convolved with the reference to a 
common PSF and the difference between the two 
was taken. The brightness of the OT in the dif- 
ference images were then measured by matching 
their PSF to the stars in the convoluted refer- 
ence images. The results were consistent with the 
RPHOT subtraction method within the uncertain- 
ties. We report here the RPHOT PSF-fitting re- 
sults with the 4 contaminating sources removed. 
Given the above comparison, the systematic er- 
ror from subtraction is smaller than the statistical 
errors and is not included in the reported uncer- 
tainties. 

The unfiltered thinned ROTSE-III CCDs have 
a peak sensitivity in the R band wavelength range. 
We thus calibrate the zero point magnitudes us- 
ing the median offset from the USNO B1.0 R band 
measurements of selected field stars. In Figure [T] 
the ROTSE-III flux densities are scaled to match 
the GROND r'-band measurements using the off- 
sets from the simultaneous temporal fit described 
in §3.4. 



2.5. GROND 

The Gamma-Ra y Burst Optical/Ne arlnfrared 
Detector (GROND, iGreiner et al"1l2008l ) mounted 
at the 2.2 m ESO/MPI telescope at LaSilla ob- 
servatory (Chile) imaged the field of GRB 081008 
simultaneously in g'r'i'z'JHK starting 13.6 ks af- 
ter To, after local sunset under clear sky condi- 
tions. A total of 28 g'r'i'z' images with integra- 
tion times of 66 s, 115 s and 375 s were obtained 
during the first night post burst. At the same 
time, the NIR channels were operated with a con- 
stant exposure of 10 s. In addition, the field was 
observed with GROND at day 2, 5 and 10 after 
the trigger. Data reduction and relative photom- 
etry was ca rried out using standard IRAF tasks 
dTodvl 1993h . s imilar to the procedure outlined in 
Kriihler et al.l (|2009h . To exclude a possible con- 
tamination of the measured afterglow brightness 
with light from the host 2.2" away, the results from 
the PSF photometry were checked against mea- 
surements taken using small apertures (0.3 and 
0.5xFWHM) where an aperture correction was 
applied. Image subtraction, using the frame taken 
10 days after the burst as reference, further con- 
firmed that the contamination is negligible. Abso- 
lute photometry was measured with respect to the 
primary Sloan standards SA107- 1006 and SA112- 
805 in g'r'i'z' ([Smith et al.ll2002l ). observed under 
photo metric conditions and 2MASS field stars in 
JHK (ISkrutskie et al.ll2006h . 

The position of the afterglow is determined 
as R.A. = 18 /l 39 m 49 s .88 Deck = -57°25'52".8 
(J2000), with a 90% uncertainty of 0.3 arcsec. We 
adopt this as the best burst localization. Along 
the line of sight to this p osition, the Galacti c 
extinction is E(B-V)=0.095 (|Schlegel et al.lll998l) . 
and the Galactic colu mn density is 7.1 x 10 20 cm~ 2 
(|Kalberla et al.ll2005h . 



3. Spectral and Temporal analysis 
3.1. Prompt Spectral Evolution 

We first examine the spectral evolution of the 
burst. Available BAT and/or XRT data are 
extracted in three time intervals, To+55 s to 
T +85 s, To+161 s to T +195 s and T +245 s 
to To+270 s (as marked in Figure [1]), centered 
on the two broad peaks detected by BAT and 
the one X-ray peak shortly afterward. Because 
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of the strong hardness ratio evolution, we limited 
the width of the chosen time intervals. The spec- 
tra were modeled with either a power-law or a 
power-law with exponential cutoff. When X-ray 
data were available, we included in the model a 
fixed Galactic absorption component and an ad- 
ditional variable system at the host redshift. The 
fit parameters for the best models are listed in 
Table Q] We tried modeling the spectra including 
a hig h energy power-la w tail using a Band func- 



tion (jBand et al.l 119931 ). but the power-law index 



above the cut-off energy is not well constrained 
by our data. The best fit spectral models at each 
epoch are plotted (with labeled offset for clarity) 
in Figure [2] as solid lines. 




0.001 0.010 



Fig. 2. — Spectral models in the selected time in- 
tervals. Interval 1, 2 and 3 (offset for clarity) cor- 
respond to the first 3 prompt peaks. Interval 4 
corresponds to the bright X-ray flare seen by the 
XRT and 5 is the X-ray afterglow. Optical fluxes 
are corrected for Galactic extinction, but are not 
included in the model fit in intervals 1 to 4, be- 
cause the optical fluxes are likely dominated by 
the external forward shock and thus not related 
to the high energy emission. We include these 
points as upper limits to constrain the internal 
emission component in the optical. The dashed 
lines, connecting the data points in the X-ray and 
the optical regions in the same interval, do not 
represent the actual spectral shape. During inter- 
val 5, the emission is entirely from the afterglow 
and the spectrum is plotted in dashed-dot line. 

In both interval 2 and 3, the high energy spectra 
are well described by a power-law with exponential 



cutoff model where the cutoff energy drops from 
above 60 keV in interval 2, to 15 keV in interval 3, 
and the photon index softens. The measured local 
column density (L06(-0.27,+0.32) x 10 22 cm~ 2 and 
1.01(-0.26,+0.29)xl0 22 cm~ 2 ) shows good agree- 
ment between intervals. During interval 1, the 
BAT data are fit by a power-law model with a 
relatively hard photon index (r ^1.33). Given 
the general trend of softening, it is likely that the 
cutoff energy at this early epoch is close to or 
above 150 keV and thus not well constrained by 
the BAT energy coverage. ISakamoto et al. (2009) 
has demonstrated that a reasonable estimate of 
the peak energy (in vf„ space) can be achieved 
using the power-law photon index in the BAT en- 
ergy range if the latter is between 1.3 and 2.3. Our 
value of 1.33 for the 1st interval is just inside this 
range and corresponds to a peak energy of 196(- 
47, +587) keV, assuming a cutoff power-law model 
as the actual spectral form. For the later two inter- 
vals, the peak energies can be explicitly calculated 
as (1 — (3)E cut , which are ~47 keV and ~6 keV 
respectively. These numbers ag ree with the es- 



timate s following the method in ISakamoto et al 



(l2009h . although the BAT- alone power-law pho- 
ton index in the 3rd interval is not well constrained 
and exceeds 2.3. 

During the 2nd and 3rd intervals, optical data 
are also available, but they are not included in the 
spectral fit for two reasons. First, emission from 
a different origin (e.g. external forward shock) is 
likely dominating the optical data (see §3.3 and 
§3.4). Second, only data in the white filter (or 
unfiltered) is available, thus any optical deficit be- 
low the extrapolation from high energy observa- 
tions can be compensated by an arbitrary value of 
local extinction, while any optical excess is hard 
to accommodate within the chosen models. In- 
stead, we plot the white band fluxes (corrected 
for Galactic extinction and plotted with the same 
offsets as the corresponding high energy spectra) 
in Figure [2] and connect them directly to their 
high energy counterparts with dashed lines, just 
to show how they compare. The optical detections 
are considered as upper limits to constrain the in- 
ternal emission component. If internal emission 
has contributed to at least part of the optical ob- 
servations, they do not exceed the extrapolations 
from XRT and BAT detections at any time. 
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Table 1 

Best-fit spectral parameters to the burst and the flares 



Interval 


OO 


Host nH 
(10 22 cm- 2 ) 


|8 


E cut 
(keV) 


X 2 /dof 


Instrument 


1 

2 
3 
4 


55 - 85 
161 - 195 
240 - 275 
345 - 385 


1 06+ ' 32 
1,uo - 0.27 

1 oi +0 - 29 

i,ui - 0.26 
, 77 +0.38 
x - ' ' -0.32 


0.33 ±0.07 
26+ - 04 

58+ - 12 
q r-,+0.13 
lol -0.12 


63.2 + * 3 8 3 

W.9±|;| 


58/56 
201/184 
200/194 

97/75 


BAT 
BAT & XRT 
BAT & XRT 

XRT 



Note. — All errors quoted correspond to 90% confidence intervals. 



3.2. X-ray Flare 

A bright single X-ray flare is superimposed on 
the rapid decay phase after the second broad X- 
ray peak. It has a similar hardness ratio evolu- 
tion as the earlier peaks. It also shows a fast-rise- 
exponential-decay (FRED) profile that is typical 
of a GRB pulse. Its rapid rise and decay supports 
an identification with extended central engine ac- 
tivity. We thus analyze the flare within the same 
frame as the prompt emission. 

We extract the X-ray spectrum from To +345 s 
to T +385 s (noted as interval 4), centering on 
the flare. UVOT v band and ROTSE-III detec- 
tions are available during this period, but there is 
no sign of a flare in the optical indicating that the 
emission is dominated by the afterglow. So we fit 
the XRT data alone. The spectrum is reasonably 
fit by a single power-law model with a photon in- 
dex of 2.5 (see Table [lj, which is much softer than 
the prompt emission and agrees with the softening 
trend. The inferred column density is somewhat 
higher than those from the earlier times, proba- 
bly caused by curvature in the real spectral shape. 
This would suggest that E pea k is near or within the 
XRT energy range at this time. 

3.3. Early Optical Afterglow 

The ROTSE-III observations show an ini- 
tially rising optical transient with two peaks at 
T + -135 s and T + ~170 s (see Figure Q] and 
[3]). The later one is almost simultaneous with the 
BAT and XRT peaks and may be attributed to a 
related prompt optical flare. The first one, how- 
ever, does not appear to be correlated with the 
high energy emission. As shown in §3.4, the over- 
all optical lightcurve is well-fit by one dominating 



component. We therefore interpret this first peak 
in terms of external shock emission. A detailed 
discussion of this interpretation will be presented 
in 84.1. 



BAT Count x 800 
XRT Count 
UVOT White x 2 
ROTSE-III ->r' 



i + 



+ T 



. +f+ , . 



1 ? + 



+ + t 



Fig. 3. — Contemporaneous BAT, XRT and opti- 
cal (ROTSE-III and UVOT white) observations. 
Time bins for the BAT and XRT data are de- 
termined by Bayesian block analysis (see text in 
§3.3). The high energy data points are connected 
by dotted lines as guidance to the eye. The UVOT 
event data is binned in 10s intervals. The smooth 
rise and decay model is plotted over the optical to 
show that no significant deviation is detected. 

Figure |3] shows a close look at the BAT, XRT 
and optical lightcurves before To+345 s. The BAT 
and XRT data are binned in intervals deter mined 
by Bayesian block analysis (|Scargld 119981 ). For 
BAT, this analysis is carried out by battblock on 
the 1 s-binned lightcurve data. For XRT, only 
WT data are used and the background (~1 ct/s) 
is negligible compared to the bright signal. We 
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Table 2 

Afterglow temporal best-fit parameters 



Energy Range 


Rise time 


Start of Decay 


ai 


Tbroak 


OC2 




t a (e) 


T„(s) 




w 




X-ray 


130* 


1800** 


0.96 ± 0.10 


15870 ± 8623 


1.78 ± 0.25 


Optical/IR 


130* 


136 ± 18 


0.96 ± 0.03 


7834 ± 2160 


1.313 ± 0.048 


X-ray & Optical/IR. *** 






0.81 ± 0.07 


15139 ± 8150 


1.71 ± 0.13 



* Fixed parameters, to a value estimated from the optical fit. 

**The plateau is hidden by the prompt emission, thus no error is estimated. 

***Only data after To + 750s are fit by a smoothly joined broken power-law model, with a smooth factor, s=l. 



Note. — All errors quoted correspond to 90% confidence intervals. 



thus adapted the algorithm used in battblock for 
event files and consider only the change points re- 
sulting in time bins longer than half a second to 
ensure comparable signal to noise ratios. The high 
energy lightcurves show significant variations on a 
similar time-scale; while no convincing short-time- 
scale variability is detected in optical. This is not 
surprising though, as the optical curve connects 
smoothly into a monotonic decline afterward. 

3.4. Mult i- wavelength Afterglow 

We adopt the two co mponent exponen t ial to 
power-law decay model in IWillingale et al~ ( 2007 , 
Eq (2) and descriptions within) to characterize the 
temporal flux evolution in X-ray and optical. In 
this model, the first component corresponds to the 
prompt emission, and the second relates to the 
afterglow emission. All optical observations are 
sufficiently fit by one afterglow component. In the 
X-rays, the second component is not revealed until 
the end of the last flare. So the rise time (t a ) of this 
second component is fixed to 130 s according to the 
optical fit. The best fit models for the afterglow 
components arc over-plotted in Figure [T] and the 
parameters are tabulated in Table [2] 

For both optical and X-ray, a later tempo- 
ral break is required in the model. The break 
times in the two energy ranges (15870±8623 s and 
7834±2160 s) are consistent at 90% confidence, 
but the final decay (o?2 in Table [2]) in the X-ray 
is significantly faster than in the optical. We 
further attempted fitting all energy bands after 
T o +750 s simultaneously with a smoothly joined 
broken power-law model (i(t)=a[(t /tf,)~ sai + 
(t i 'if,) - ' sa2 ] 1 1 /s ) , with a small smoothness factor, 



s~l. Although the overall fit is acceptable, giving 
a reduced Chi-squared of 1.17 for 239 degrees of 
freedom (see Table [5] for the best fit parameters), 
the fit parameters are dominantly determined by 
the well-confined optical data. Some systematic 
deviations from the model, not reflected in the 
Chi-squared statistics, are noticed in the X-rays. 
Between T +7 ks and T +28 ks, the X-ray points 
are mostly above the model, and after To+30 ks, 
the X-ray data are all below the model. Our 
data are thus consistent with displaying a smooth 
achromatic temporal break, but a chromatic be- 
havior is not confidently ruled out in case of a 
subtle evolution. 

We then examine the spectral energy distri- 
butions (SEDs) in selected afterglow epochs. In 
the X-rays, no significant spectral or hardness 
ratio changes can be detected. In the optical, 
the UVOT detections before the break and the 
GROND observations after the break are poorly 
fit by one power-law spectral index plus the same 
amount of local extinction. This is caused either 
by spectral evolution or a spectral break at the red 
(blue) end of the UVOT (GROND) coverage, at 
around 2.5 eV. This latter conjecture is confirmed 
by further analysis. 

The overall shapes of the SEDs are modeled 
from optical/IR to X-ray with all available data 
before and after the break. Assuming they belong 
to the same synchrotron emission component, we 
allow either one or no spectral break within the 
energy range. Before the break, we extract an X- 
ray spectrum in the time interval from To+765 s 
to To+1778 s (first orbit) and from T +5286 s to 
Tq+7562 s (second orbit) to maximize the signal 
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Table 3 

Best-fit afterglow spectral parameters 



T-T 


Energy Range 


Host Type Host Extinction nH 


Port 




0* 


X 2 /dof 


(ks) 




E(B-V) (10 22 cm -2 ) 




(cV) 







1.5 


XRT only 






47+0-36 
u -*' -0.33 








1 06+ 014 


46/43 


1.5 


XRT & UVOT 


MW 


n 1 5+ 04 


u - DZ -0.23 


1 


.14 ± 0.05 




1.14 ± 0.05 


52/46 


1.5 


XRT & UVOT 


LMC 


n+ 02 


u.oa_ 22 


1 


.12 ± 0.04 




1.12 ± 0.04 


53/46 


1.5 


XRT & UVOT 


SMC 


0.08 ± 0.02 


54+0- 24 


1 


.10 ± 0.04 




1.10 ± 0.04 


55/46 


29.5 


XRT, UVOT & GROND 


MW 


n ,o+0 09 
0.ls_o. 07 


Q7+ 85 




n ,, + 0.30 
U-l.3_o.41 


2.4 ± 0.2 


1 12+ 011 


11/14 


1.5 & 29.5 


XRT, UVOT & GROND 


MW 


16+ 004 


0.65I 2i 




06 +0 - 17 


2 4+°- 2 


1.14 ± 0.05 


72/63 



Note. — All errors quoted correspond to 90% confidence intervals. 



to noise ratio, given that no significant spectral 
or hardness ratio changes is detected in the X- 
rays. This spectrum is scaled to the X-ray count 
rate at T +1.5 ks to minimize the interpolation 
of the UVOT data. The UVOT mean count rates 
at To+1.5 ks are estimated by fitting the nearby 
data points (within 500 s on each side) in each 
band with power-laws. The errors are estimated 
using a bootstrap approach, where each data point 
is randomly re-sampled from a Gaussian distribu- 
tion centered on the observed count rate with a 
sigma equal to the measured error. The UVOT 
count rates are converted to XSPEC compatible 
spectral files with the corresponding response ma- 
trices downloaded from the HEASARC Calibra- 
tion Database (version 20041 120vl04). 

The SED can be fit by a single power-law model 
with fixed Galactic dust and gas absorption and 
additional dust and gas absorption at the host red- 
shift (see Table . This model is supported by 
the well-matched initial temporal power-law decay 
indices {a\ ~0.96) in the optical and the X-ray. 
The relatively small amount of host extinction re- 
quired and the lack of coverage at the rest-frame 
2 175 A Galactic absorption feature do not allow 
us to distinguish between the Milky Way (MW), 
Large Magellanic Cloud (LMC) or Small Magel- 
lanic Cloud (SMC) environment. 

After the break, another X-ray spectrum is con- 
structed between To+16 ks and the last detec- 
tion at To+400 ks. This spectrum alone has a 
relatively low signal to noise. We scale this X- 
ray spectrum to the latest epoch where the after- 
glow was detected in the J, H and K bands (at 
Tq+^29.5 ks) and convert the contemporaneous 



GROND detections to XSPEC compatible files. 
The UVOT u and b band rates are also interpo- 
lated to this epoch by fitting the neighboring two 
detections (about 30 ks apart) with power-laws as 
described above. 

A single power law model does not produce an 
acceptable fit to the overall SED, regardless of 
the choice of local dust and gas absorption. At 
least one spectral break is required between X- 
ray and optical. Our best fit model (see Table \5§ 
constrains a break at ^2.4 eV, in the UVOT v- 
band. Such a low break frequency is not well cov- 
ered by the earlier UVOT data, and the two SEDs 
thus have consistent shapes. A MW extinction 
curve is used because it provides a best fit and 
is marginally preferred during the earlier epoch. 
Nevertheless, we do not consider this as a conclu- 
sive property of the environment because of the 
possible uncertainty in our selected intrinsic spec- 
tral shape, e.g. the smoothness of the break is not 
considered in our model. 

Figure [4] shows a simultaneous fit of the two 
SEDs with a broken power-law model with Galac- 
tic dust and gas absorption and dust and gas ab- 
sorption at the host redshift. We tie the dust (fol- 
lowing a MW extinction law) and gas properties 
for a consistent environment throughout the af- 
terglow phase. The fit is fairly good except in the 
uvwl band. A plausible explanation for the under- 
estimate in uvwl is that our model applies a sharp 
cutoff above the Lyman-break which might not be 
true. The best fit model parameters are consistent 
within errors with the estimates obtained in nar- 
rower energy bands and smaller time coverages as 
described above. Also, the column density derived 
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is consistent with the estimate during prompt in- 
terval 2 and 3 (see Table [T]) . 



could have increased over time from below optical 
to the observed 2.4 eV. 




Fig. 4.— Combined XRT, UVOT and GROND 
spectral energy analysis during the afterglow 
(epoch 5 and 6). The local extinction and absorp- 
tion properties are tied. The fit parameters are 
listed in Table |3l The absorbed models are over- 
plotted in the upper panel. Dotted lines show the 
intrinsic power-law and broken-power-law compo- 
nents. The location of a spectral break is indicated 
by an upward arrow. The residual ratios are plot- 
ted in the bottom panels. 

In summary, the observed optical and X-ray 
afterglow emission started off following similar 
temporal decay rates (a\ ^0.96) and displayed 
smooth temporal breaks at similar times (at 
~15 ks). While the two SEDs before and after the 
break were characterized intrinsically by a power- 
law model and a broken power-law model respec- 
tively, they had consistent shapes because the low 
break frequency was not well covered during the 
first epoch. Besides, the optical/IR observations 
by GROND between T +13 ks and T +30 ks do 
not show evidence of a moving spectral break. The 
overall afterglow behavior was consistent with be- 
ing achromatic. Nonetheless, there is some hint 
of a slightly steeper decay in the X-rays after the 
break. Because the afterglow was not observed in 
IR and near-IR at early times, probable evolution 
of the spectral break is not constrained, e.g. it 



4. Discussion 

4.1. Early Optical Afterglow 

4-. 1.1. Reverse or Forward Shock 

For the initial optical observations, fixing the 
beginning of the external shock to the burst time 
To yields a temporal rise index of a = —5.0 ± 0.6. 
This is too rapid to be considered as the pre- 
deceleration phase of the forward shock (i 2 or t 3 
in ISM, slower than t 1 / 2 in Wind medium, see 
Panaitescu fc Vestrandl200"8h . but agrees well with 
the prediction for a reverse shock rising in a thin 
shell case when the burst duration is less than the 
decel eration time ( Kobayashi 2000t Zhang et al 
20031 ) . Back-extrapolation from late time after- 



glow observation shows that the reverse shock and 
forward shock must have comparable peak lumi- 
nosities. The analysis of the later afterglow emis- 
sion suggests the forward shock characteristic syn- 
chrotron frequency (y m ) may have passed below 
the optical by T o +200 s. If the micro-physical 
parameters are comparable in the reverse and 
forward shock regions, the reverse shock should 
have a much lower characte ristic synchrotron fre- 
quency ( Sari fc Piranl 1999) and thus have a very 
weak optical emission. The detection of reverse 
shock emission would then poi nt to a more mag- 



netized reverse shock region |"Zhang et al 



2003; 



Kumar fc Panaitescull2003t Ijin fc Fanll2007t ). 

However, there are several questions to be asked 
for this interpretation. First, is the thin shell as- 
sumption correct? Although the overall prompt 
BAT and XRT detections extend to after the be- 
ginning of the optical detections, the second epoch 
of strong burst activity emerged once the initial 
rapid rise in the optical had stopped. The tempo- 
ral and spectral similarity between the flares indi- 
cate that the second epoch of high energy emission 
was from the same emitting region as the earlier 
bursts of emission, i.e. closer to the central en- 
gine than external medium. This later part of the 
ejecta would then run into the external medium 
after the initial shock has been established. We 
thus can ignore the effect of this later injected 
kinetic energy when considering the early rising 
lightcurve. Suppose the bulk of kinetic energy was 
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injected into the external medium by the materi- 
als producing the first flare ending at To+100 s, 
the shock can be reasonably considered impulsive. 

This leads to the second question: whether 
the onset of the af t erglow is chosen correctly as 
To? [Quimbv et al.l (|2006l ) have pointed out that 
different choices of time origin may lead to dif- 
ferent interpretations of the early data (also see 
Tagliaferri et al. 20051 for a discus sion in the X- 
ravs). lLazzati fc Begelman ( 20061 ) explored how 
varying central engine activity would affect the 
self-similar evolution in the early afterglow and its 
lightcurve. Their simulation is limited to a power- 
law energy release, but gives a good schematic pic- 
ture. They found that if the burst energy release is 
concentrated toward the end of the burst duration, 
the self-similar evolution would deviate from one 
simple power-law behavior. For GRB 081008, the 
emission at T +65 s is about one order of mag- 
nitude higher than at To, it is therefore proper 
to choose a delayed time origin in a power-law 
fit. Such a small time shift has a negligible ef- 
fect on the temporal index inferred for the later 
smooth decay. But if we choose the reference time 
to be To+45 s, the initial rise phase would have 
a temporal index of —3.1 ± 0.4, fully consistent 
with the pre-deceleration of the forward shock. 
For any choice of reference between To +25 s and 
To +62 s (including the beginning of the first big 
pulse detected by BAT), the estimated tempo- 
ral power-law index is consistent with -3 within 
2 sigma. Given the earlier discussion about the 
reverse shock and the fact that the lightcurve is 
well fit by one smooth rise and fall component, we 
prefer the external forward shock interpretation 
for this rising lightcurve observed by ROTSE-III. 

Assuming the first emission peak in the optical 
marks the onset of the forward shock deceleration, 
we can apply the adjusted rise time (To+135 s - 
T +45 s) of 90 s and the energy released before 
Tp+1 35 s (see §4.1.2) into Eq (1) in lMolinari et al ' 



(|2007l ) to obtain a Lorentz factor of Tdoc —225 
at the deceleration radius a nd an initial Lo rentz 
factor T of -450 (~ 2r doc ; iMeszarodExih . We 
have assumed a particle density n=l cm -3 and a 
radiative efficiency /i=0.2. The deceleration radius 
is then estimated as — 2r 2 ci— 3 x 10 17 cm. 



4.1.2. Refreshed Shock? 

After the forward shock start decelerating, the 
passing of the synchrotron frequency and/or cool- 
ing frequency across the optical band may produce 
temporal breaks. But a second peak, as observed 
for GRB 081008, is hard to explain in a model 
involving only one external forward shock compo- 
nent. The temporal coincidence with a high en- 
ergy flare suggests possible contribution from in- 
ternal emission. This scenario is hard to rule out 
given the slight optical excess (—4 sigma) above 
the afterglow model seen in Figure [3] However, 
the lightcurve also shows a probable dip below the 
smooth, one component model, and the decay af- 
ter the second peak joins smoothly into the later 
phase. We thus explore other interpretations for 
the double peaked lightcurve. 

In a reasonable contemplation in the frame of 
the fireball model, the shells of material produc- 
ing the late burst activity would reach the external 
medium and catch up with the early decelerated 
shells at some moment and re-energize the forward 
shock. This later scenario results in an irregular 
"step" structure, as additional energy is injected 
over a time period comparable or longer than 
the second epoch of burst activity. It has been 
noted that refreshed shocks from slow or late shells 
may alter the shape of the afterglow lightcu rves 
(|Rees fc Meszaroslll998HKumar fc Piranll2000l ) . It 
is, however, hard to directly relate such structure 
with the internal emission as they are produced 
by different mechanisms at different regions. Here, 
we provide a simplest possible test on this specu- 
lation. 

If a second batch of kinetic energy is injected 
into the forward shock, a new self-similar solution 
would be established sometime after the cessation 
of all burst activity. Indeed, we didn't observe any 
apparent deviation from a power-law decay after 
To+200 s. We try to estimate the synchrotron 
emission before and after the refreshed shocks. 
Without the additional energy injection, we can 
assume the afterglow would decay monotonically 
after the first peak at the same rate as observed 
in the later afterglow. At a time after To+200 s, 
the observed flux is 50-60% higher than this ex- 
trapolation, including uncertainties from the flux 
measurement at the first peak and the normaliza- 
tion of the later decaying lightcurve. 
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We then estimate the energy release in the 
prompt emission before and after To+135 s (the 
first optical peak). We estimate the fluence be- 
tween 0.1 keV and 10 MeV in the rest frame. In 
the observer's frame, the 0.1 keV lower bound fits 
nicely within the XRT sensitivity range. Before 
Tq+135 s, only BAT observations are available 
and the average spectrum can be fit by a power- 
law with exponential cutoff model. The cutoff en- 
ergy at 157 keV, although above the BAT energy 
threshold, agrees well with calcu lation using the 



power -law photon index following Saka moto et al 
( 20091 ). We estimate the fluence during this pe- 
riod to be ~ 4.7 x 10~ 6 erg cm" 2 . The hard 
X-ray emission dropped below the BAT detec- 
tion threshold at To+^265 s, but it is likely that 
the central engine activity continued with softer 
emission detected mainly in soft X-rays. Be- 
tween To+155 s and T o +210 s, simultaneous XRT 
and BAT detections constrain the spectral shape 
well, and using a power-law model with an ex- 
ponential cutoff, we estimate a total fluence of 
r~> 2.5 x 10~ 6 erg cm~ 2 . After this, the peak en- 
ergy shifts further toward lower energies and the 
fluence in the next 100s is only one fourth of the 
previous interval (with cutoff energy at ~10 keV). 
The additional energy release between T +135 s 
and To+155 s is even smaller. From T +135 s to 
To+^310 s, we estimate a total prompt fluence of 
- 3.4 x 10~ 6 erg cm- 2 . This is about 70% of the 
total energy released before To+135 s, and consis- 
tent with the flux increase calculated within the 
large uncertainties. Note that the overall energy 
output is about 20% higher than the estimate in 
§2.1 using an average gamma-ray spectrum and 
the BAT detected burst duration. 

The X-ray flare at To+~400 s is likely gener- 
ated by the same internal mechanism. We can 
estimate its emission in the same energy range as 
above. The spectrum is so soft (with photon in- 
dex greater than 2) that we can directly integrate 
it in the high energy range. No extrapolation be- 
low the XRT energy is necessary with our choice 
of energy range, so we do not have to worry about 
a spectral break below 0.3 keV. The "bolometric" 
fluence calculated is two orders of magnitude lower 
than the previous prompt emission. Therefore, if 
the contribution to the forward shock is propor- 
tional to the prompt luminosity, the change in the 
forward shock emission is comparable to the error 



bars in the optical measurements and too small to 
be noticed. 

A caveat for this calculation is the assumption 
that the relative radiative efficiency in the internal 
and external processes is constant although the re- 
lationship between the two is not straightforward. 
During the burst, the peak emission frequency 
down-shifted from 7-ray to X-ray. The X-ray emis- 
sion could be from slower or delayed ejecta, or al- 
ternatively from collisions with a smaller Lorentz 
factor dispersion and thus lower efficiency in con- 
verting the kinetic energy. At the forward shock 
front, the conversion efficiency may depend on the 
condition of the previously shocked medium and 
therefore requires a detailed simulation to deter- 
mine. Another overlooked factor is the evolution 
of the synchrotron spectral shape, as the flux is 
estimated in a narrow optical band. 
( 2006h 



Vestrand et al 



showed that the early 
optical lightcurve of GRB 050820A can be mod- 
eled with a prompt component tracing the 7-ray 
emission and two afterglow components. The af- 
terglow components, correspondingly produced by 
the precursor and the main event ejecta, have con- 
sistent ratios between the prompt fluence and af- 
terglow amplitude. This agrees with our finding 
although their first afterglow component is faint 
and not well-sampled for a detailed study. A sim- 
ilar double peak ed prompt lightcur ve is observed 
for GRB 061007 (|Rvkoff et al.ll2009h . The overlap- 



ping high energy pulses indicate continuous energy 
input into the external shock and thus modeling 
of the early lightcurve in a thick shell case may be 
more appropriate. For that event, the initial very 
rapid rise of the optical emission is also problem- 
atic in the current forward shock model. 

4.2. Late Afterglow Evolution 

Before the temporal break and above —2.4 eV, 
optical and X-ray appear on the same segment 



of 



1 sp 
intE 



ectrum, we 



a synchrotron 
method described 
all scenarios consistent with 



Racusin et al 



thus use the 
(|2009l) to find 
the observed tem- 
poral and spectral indices. The temporal decay 
(a=0.96±0.03) marginally agrees with the value 
expected ((3/3-l)/2~1.2, for /3 -1.1) for above the 
cooling and characteristic synchrotron frequencies, 
regardless of slow or fast cooling regime and the 
external medium density profile. In such a sce- 
nario, both the cooling frequency and the charac- 
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teristic frequency are below or close to the optical 
by To+200 s. 

For all the other closure relations considered, 
the temporal decay is shallower than predicted if 
no energy injection is prese nt. In terms of the 
canonical X-ray lightcurve (Nousek et al.l [2006; 
Zhang et al.ll2006UO'Brien et alJl2006h . this decay 
index is within the range observed for the shal- 
low decay phase (phase II). This phase, detected 
in a significant fraction of the bursts, is generally 
explained as due to prolonged energy injection. 
However, the source of energy is unclear as well 
as the mechanism for fine tuning the input energy 
into the forward shock smoothly over such long 
period of time. 



We apply the evolution constraints in lRacusin et al 



( 20091 ) to the X-ray data and find the cause of the 
break is consistent with either a jet geometry or 
the cessation of energy injection. For both types 
of break, achromatic behavior extending to the 
optical range is expected, consistent with our ob- 
servations. 

After the temporal break, our data tentatively 
suggest a spectral break around 2.4 eV, although 
it is hard to constrain the exact spectral shape 
and the environment properties at the same time. 
Such a break is hard to be accommodated in the 
standard fireball model. A synchrotron spectrum 
is characterized by three possible break frequen- 
cies, a cooling frequency (v c ), a characteristic fre- 
que ncy (v m ) and a self- absorption frequency (v a ) 
(e.g. lSari fc Piranlll999t[Granot fc Sarill2002h . All 
three frequencies should evolve over time. The 
observed break at 2.4 eV is too high to be a self- 
absorption frequency Given a spectral index be- 
low the break (-0.13{-0.41,+0.30}) consistent with 
-1/3, the shape of the spectrum can be expected 
i n a slowing cooling regime (y m < v c ) around v m 
dGranot fc Sarill2002h . However, v m should always 
decease over time as the forward shock decelerates, 
unless the shock is reenergized by significant late- 
time injection, which does not seem to be the case 
here. We do not observe the expected rise or flat 
evolution in IR (below the break) and do not de- 
tect a higher break frequency in our earlier SED. 
Finally, our data do not rule out an increasing 
break frequency, as expected for a cool ing break 
in a wind medium ( Chevalier fc Li 2000j), but the 
spectral index difference across the break (~1) is 
much larger than the predicted amount of 0.5. 



If some chromatic effect exists around the tem- 
poral break, the passage of a cooling frequency 
in a wind medium produces the correct trend of a 
break moving from optical to X-ray. The expected 
temporal index change of 0.25 is similar to the ob- 
served value in the optical (~0.3) but not enough 
to account for the amount of change detected in 
the X-rays (~0.8). In any case, the temporal break 
at T +15 ks is not likely to be caused by the move- 
ment of the observed spectral break, because there 
is no evidence that the break frequency has yet 
reached the X-ray regime by T +29.5ks, as would 
be expected in such a case. 

In short, the observed temporal break is con- 
sistent with being a jet break or the cessation 
of energy injection. Before the break, the multi- 
wavelength observations are best interpreted by 
the cooling and characteristic synchrotron fre- 
quency being below the optical. No satisfactory 
explanation is found for the later spectral shape 
and its probable evolution. It is viable to assume 
that the spectral break has increased over time, 
as the cooling frequency in a wind medium, but 
such evolution is not likely to be responsible for 
the observed temporal break. 

4.3. Burst Energetics 

If the temporal break at To+^15 ks is inter- 
preted as a jet break, the jet opening angle can 
be estimated as 2.1° or 2.4° for ISM dSari et al' 



1999) or a wind medium (IBloom et al 



20031) re 



spectively. In the calculations, we have assumed a 
particle density n=lcm~ 3 , a characteristic wind 
density A*=l, a radiative efficiency ^=0.2 and 
a kinetic energy Ef. = Ei so //i. The collimation 
corrected total gamma-ray energy release (E 7 ) is 
hence 4.2 x 10 49 erg and 5.5 x 10 49 erg for ISM and 
wind medium respectively. This number would be 
slightly larger if early X-ray detections are also 
considered (see §4.1.2). 

For a jet break, the change of the t emporal in- 
dex ( ~0.8) favors an ISM- like medium ( Panaitescul 
20051 ) . Both the inferred opening angle and the 



collimated energy output are on the small side 
but well wi thin the distribution for Swift detecte d 
GRBs fe.g. lLiang et alJl2008llRacusin et al.ll2009l) . 
The nature of the break is, however, uncertain 
given the previous discussions on a probable sec- 
ondary chromatic effect. A later jet break, occur- 
ring close to or shortly after the afterglow dropped 
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below the detection threshold of all our instru- 
ments (^400 ks), would also give a sensible es- 
timate of the opening angle and E 7 . 

5. Conclusions 

With energetics typical among Swift detected 
events, GRB 081008 happens to last long enough 
to have overlapping phases with the early follow- 
up observations in X-ray and optical. It is well be- 
haved in the sense that all early canonical phases 
in X-ray and optical were detected. A key feature 
of the high energy emission is the hardness ratio 
evolution in accordance with the flux fluctuation 
throughout the burst and the X-ray flares. The 
temporal continuity and spectral softening suggest 
that early X-ray observations originate from ex- 
tended internal activities, but either with a grad- 
ually decreasing efficiency or a declining energy 
output from the central engine. 

In the optical, the detections seem to be dom- 
inated by afterglow emission at all times. An in- 
triguing question to ask is what happens to the 
material producing the late burst emission. In 
the frame of the fireball model, they will run 
into the ambient medium and refresh the exter- 
nal shock. Whether this signature is observed for 
GRB 081008 is uncertain. If such early refreshed 
shocks can be observed, their relative strength to 
the prompt emission provide important clues to 
the burst mechanism and production of the after- 
glow. Refreshed shocks have been used to explain 
late-time lightcurve bum ps ( Granot et al.l 120031 : 



Johanncs son et al. 200G). In those cases, slow 



shells catch up with the forward shock front hours 
to days after the initial burst. The existence and 
properties of the slow shells are only inferred from 
the lightcurve shapes. Observations directly link- 
ing the ejecta material and refreshed shocks are 
achievable for long bursts with discrete episodes 
of activities or for events with very bright X-ray 
flares if they are produced in a similar way as the 
prompt burst emission. The latter cases may pro- 
duce cleaner chasers of the refreshed shock with 
less contamin ation from internal o ptical emission 
(howe ver, see Kriihler et al. 20091 ). Falcone et al 



( 20061 ) have attempted to relate the subsequent 
bumps in the X-ray lightcurve with the bright X- 
ray flare in GRB 050502B. No solid evidence of 
refreshed shocks associated with flares has been 



observed in the optical so far, probably due to 
limited time coverage of the data or incomplete 
understanding of the modification to the afterglow 
lightcurve. 

The importance of well-sampled multi- wavelength 
data has been demonstrated. In general, the ob- 
served afterglow emission is consistent with dis- 
playing an achromatic temporal break. Such a 
break can be interpreted as due to a collimated 
outflow (jet) or the end of energy injection. How- 
ever, there is a hint of a steeper decay in the 
X-rays after the break. While for GRB 081008, 
the deviation is rather subtle, similar steeper de- 
cay in the X-rays than in the opti cal has been 
observed for some other events (e.g. iRvkoff et al.l 



2006tlPage et~aHl2009l) . Using data from GROND, 
UVOT and XRT, we are able to model the SED 
of the afterglow from IR to X-ray and constrain a 
probable spectral break within the optical regime. 
As many other well-studied bursts, the overall 
temporal and spectral behavior of the afterglow 
is hard to explain in the frame of the standard 
fireball model. 
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